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و۱۲ ما۸ موی رمع مصلو۲] 26۵ وزمع۱ ممل) ۳۵۵60۵ امد عمتام ۱۷۲00 
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ممامدنسد 2« ,۲۱ ,ازمعطمی1 ۲ .1۷۲ زنمطووم۸ 51۰ 


0ج اتود گه لمح رصتهعصتعصظ مه وعمممزهو ار50 0۶ )ممحص 1۳۵۵2۲ ,۳۶۵۲6950۲ ۸99001216 220 ۲۵۶6990۲ -3 20 1 
۰ ر1121 ,۸02011 رتللطاهت۸ طومطمممم ۵۶ اهتنا ردعمتتامفم۴ ۱۱26۲21 

(26.15 :۱۱۵۵ ۵) زعقطاع وقطه :تم عمطانم مصتطممومترمن -۶) 

۶ الوه۷تصت۲] روعمامومک لهدل مه امه له لنمه۲ راصمصصندمونا مصلهعصتمص تماه ۷۷ ,۲۳۵۲6990۲ )صماوزوو۸ -2 
20 ,۸۲02011 ,۸۳020111 طعمطعهطه]۷( 


:6 نها ان 10 ۲۱۵۷ 12004 360۷۵ 
5011 ۱۷۲۵0۵1102 2024(۰) ۲۱۰ ,0م21قط)آنه حاقطقطای ی ر. متصهطعهکا جیامموصهوه۲۱ بیط رتتفطعفخظر 0 ۵ 
2676 270 ۱6۱۷0۲1۲ متام لملمتتا نم رجمتووم‌نعع۲ عملدنا ممصه) فلع ممتاقتاعص۵ه 19-03-2054 : 


:۵0 ۸ 
اعتاعص طاز۳ صعتفم۴ ط1) ,269-283 ب(38)2 رات ۵۳۸0 ۲۷۷۵۸۵ ۵۴ م۵ .عصتصهع2۳0۵ 19-03-2034 ۳ 
85 ۷ ۰0۲2/۸ 05://001ا۸ظ .(امعادتاج ۱ 
9( 


۲ 1386 )5و2 [0ز عط) 0۶ ممصماولوع لهعنصهطممصه فطع و«مطاو لصو عط ۵ (کزظ) معصه)ولوع۲ ممتاقتاعجمم مط 1 
ب60611ع م1 .فصمتاهتهون مات هه ظ)2۲09۷ )۲۵0 مصمتافصتصهع 9660 01 معا طا کصهاز0 م1۳ ور 1 رعما0۲0 ال عم آهعنص۵ع 0۴ 
۰ )مهم تصوژو 0۵ ۷۷۱11 5011 مطا ما فامم ۵ ممتوصحمه هه طا7«مع معط ب۵ظ۱۷]۳ 2.5 ۵۵۵605 8011 2 0۶ مباله۷ ۲ 2 ]1 
مدا 15 1 ,۲۳۵۵۲۵۲01۵ .قرمتم لو)صمحطرتا تافص معا فییل اقع) لادم صه فیامتع120 ج مقله وا کاظ ۵۶ )ممرمتنافقمجظط عع17 مط ]1 
مصمع 0ج لاله ۵60۵ همه لمتمتاه ,را ممتووم‌نوع تومصتا مامت‌لنهد مه تاو وامعمقط عصمععلن ۵۶ عفن فص 
6 .و6)وهامه هط 5011 ]1۵0۷-۵۵5 مضه فاطاتووععع2 بانقده طوبمعط کاظ عامصصتاوع ما (رظنن) و۳0۵ مملففعتمرم 
7 مطا مگ کاظ مصتامصتاوه ۶۵۲ داعم ی مه لاله یکلا صنهاطاه م۵ (1) :۷۵۲۵ طمنمعوعد ونطة ۵۶ وع00۷عز00 
0 (02 ,0۳0۷1۵68 ۸02011 ۵۶ جمنعع مملطاعم0صه۴۲ ۶ وصق1 ۷260 لته جع معصه: راوعع10 ما وماطهتنه۷ 011و 2608581016 
ممتامصتصصتعاع 0۴ اممتمتللهمه معط عصتویه کاظ ارمو عصتامصصتاوع صا واع0‌مصه 0مصمتاصمصصعتم2 فطا ۵۶ «مه‌تتاهع2 معط ععهمعمم 
تمعن رگا )همه عرکتآمابک مدآ 0ص رظا) عمته صفعصه ررظ۲۱۷/6) عم عتمناوه مصقعصه اموم رح) 


1۷۲26۳۵216 20 ۵5 


0600۳ 011و جج 0-10 حصمظ ملق «اامتامصماوره اوه مهس (ر80 ع ص) ومامرصهه عاونا صرح 1156۵60 

۶ 109عع۳ وماطعم0صه۳ ۵1 وقصعا (37 < 6 ۵۵0۷۵260 200 (23 < 08 ۲2820 ر(20 ع 68 ]10169 صا معصقافتل و 50 2۵2۲1۷ 
02۰6 ۱۵۲ص( 57 33 48۶ ۵ 457 32 48٩‏ .1082 رل 257 24 38٩ 24 10۳ 1۵0 38٩‏ .126) صق1 ,م0۲0۷ ۲02011 
,(00) ط0هاده متصهعع0 ,(رباظ) «اتدصعل مامتانهم 20 (ناظ) اه ۲ رم رون 2 هه رالد رصق ۵۶ فاصعنجمع 
0 ۱/۵۲۵ (نز]0(۷) عون موعمه متامصممع 20 (۲ ۷۲۷۷ مامصصهنل اطعله۳۷ صقمجظر رر۳۷۷) امماجمن تمه 11610 
مهن متتامموهمع صهع]۱۷ بعهل هه 0ص ناظ وتو 0ماقمانملهه و72 (امناظ) اتفصعل 0۵11 ۷۵ 20ام .ورماه12001 مط ط1 
6 .006۵01282608 0127 0ص2 ]1رد رصهه رها 0مانام‌رمی مم۷ ومآمتانهم 5۵11 ۵۲ (یم) صمت00712 270صهاو متتاعصط0عع 20 (م) 
۰ د 24 ([۵06ظ وملمصه) ععامعمههمعصوم عصمی تون حالو طاً 2628۲۵0 مه 5011 عطع ۶ (کزظ) ممصم)ولدع۲ ممتاحتامصهم 
2 99و م1 .۵0616 ۵۶ عصتاوعا 10۲ 02162 20 40ص عصتصتهتا 10۶ 0242 60 و2 وعکرع5 ۲۵ صاً 0۱۷1060 ۱/۵۲۵ «اصطرمل‌حنه؟ 12۸۵ 
تقمصنا مامتالناحظ 0۵۷۵ ما ۱۵۵0 ۷۷۵۲۵ ۵1]۷۵۲۵ 4:0۵ ۰۲۵0۵18 2٩0۲0‏ 0686 0ج هشن ۱۷۲۸1 ,مطامجط ۹060۳۷186 ۷16 ٩01100۷2۲6‏ 
۰ ,۵00615 (ظظ۳)) عصتصصهتو ۵۲0 ۵۲۵8910 عصعع مه (اآانض) ار۱۵00 تمتیامه لمل تاه ,رک اا) صمنووعتوع۲ 
تمه 10مصع 1و کصم‌عصه] مط مه کاع0اعه ممامم۵۵۲0 (123/6۲ عبط ما فصمعتامط 6 20 5 ر2) 3766-12۷6۲ 10۲۷۷2۲0 1660 خر 
طیان اوه ۵ ممتمم06۷۵1 06۲0۲۵ معومهه ۷۷۵۲۵ وتعامصصه‌تهم لمصتامن 01 ]56 ۸ .عصتام0مه آحآاخ معط ۶0۶ ۱960 ۷۷/۵۲۵ 6)108] 
۲ هه اه مطا رها 9۵۱60۲۵0 ۷۷۵۲۵ ممتامعصن وصتاصنا همه 5126 240عظ روعصعع 204 0۵5ومصصمعطن ۵۶ تعجامصناه فظ 1 ,2۵061 


فصمصصصومت و0۵2۲ تملمنا 0مباطاتتونل مامتایه وعععععة صهون صه_ وا فتط1. .(ق)۲مطاج 18 2024 
202021001۱ معصتتهطه رمع فاتصمهم طمتطاه ب(4۵ ۷ 00 عفصوعنا لفعمتافصعامآ 40 متاخ 
ما 601 )200۴0011۵ 21۷6 ۷۵۷۲ 26 1088 25 رقم زن حصتاتعص رصح صا مرمنامنا0 ۲۵۵۲0 قصه صمتا ماوت 23 
۰ 3۵ 20 (۲)8مطبیه لفصتعتره معط 
5 ۵ /105://001.0۲2 


۰ نشریه آب و خاک جلد ۰۳۸ شماره ۲. خرداد - تیر ۱۴۰۳ 


عتحعانا ما عصت20001 معومان ۷۵۲۵ قعماهت6و۵ عتاعصعع 0۶ م۲266 16 .راه 1۵5۵6001۷ + 20 ر8 ر3 ,30 276 «عط و2 رل‌مطامظ 
۶ امم)وهه ها امه م۷۵۲۵ ۲۴ مصتامصتهء صا فاممصر نی مد اه یکلنتا ۵۶ بمدتتهعه مط1 .وعتناه 
(6۲5 نمعئه‌تالهم عنامایکطمدا هه رظ) تمه مفعصه ررظک۵/6) مه معمنود صفعم ام ,رک صمتافمتصهاعل 

1۳ 


«عی 6( 200 165ا6ع1 


7۲ ,(13 ح< 0) 102 ب(12 ح ه) صوم۱ رصقو بر6 ح 6 صعما رهام مضه 1 ح ه) صصوم1 هآ 20ظ وتو 4160باو مفط 1 
,(66 64:79 60 13.14) 5200 0۶ ممباه۷ 1۳6 .وعووهاه مرها ر23 ع م) صوم1 اه مه را حه) چهان هلو بر14 ع م) صوم۱ رهآه 
۵4 ) داظ ,ناجهع0۵ دعقم 50.47 10 11.58) ۳۷۷ ,(۵ 7.17 ۲0 1.01) 00 ,(70 42.18 0) 2.95) 127 ,(۹6۵ 74.96 10 21.11) ]811 
۵ ۵6 .جمنوع: 0عننه عطا ۶و ملذمو عطا صذ مصمتاهتنه۷ 00مع 07760 (۱۲۳۵ 5.83 مه 1.03) ۳ 4ص (تصصم ع 1.42 ما 
عطا ما ف۲ 0365 )وم هه ( 0۵6 ارو 04۹ 6۶ ۳۷۷ طازس ۳۳ مععساعها مصمتاهآمجمه )صم‌تانصع ژد 0صبام] 
3 ,دالحتعصهمن ,۳۳ متمصنامه ما عاطافتنه غاورصدً حق وه عفن امط مه۷ وه عطا ر( 0.84 وم طانه ناه که واتمعمناام‌تاآنهر 
6 .۷۵۲1۵0168 5011 ع01م2۷71 ۲6۵0117 0ع۲باوحعجط مصرم کاظ متمصتاوه ما 0مامنا8ه0ع ۷۷۵۲۵ ۵0615جهر ری مه ۸۲۲ بک ]1۷ 
٩11] 0‏ ,)۲۱۷۷ ۷۷۵۲۵ ۲۴ عاهصصتایی م6 کماماه۷2۲ ۱162016 )عمجطر عط غقط) ۱7۵0«مو قاع0مجظ رن 0صع ۸ م۱۷ ۵۶ عتلناوع] 
هصح ,0.36 0صج ۱۲۳۵ 0.19 ,۷۲۳۴2 1.19 ,0.44 لمییوه 0مصنم‌داه مغ« متنعاتن ۲٩‏ مج ظ رظ9] ,۲۳2 ۵۶ معبلع* م1 .1ع1۱7ظ 
,۵0616 طن) جرج ۸۵۱۵۲ یکلا )عع عطاٌ 10۶ 0.03 ب۷۴2 0.13 ب۷۴ 0.91 ,0.79 ,0.92 20 ۱۷۲۳۴۵ 0.05- ,۷۲۳۴ 0.41 ,0.92 
۴ ۵60۷۵۵۵ ومتهاهزمه اصهملصع زو هه 0682076 2 و1 67 تفا ۲۵۵0۲۲۵۵ ولج فتم‌طمته6و۲۵ تمصر1۵ م1 .۲۵9۵6011۷۵1 
۷۰۰ ۷/1۳ 


صمتعصل(10) 


۲ 36 مر۱۷۵ (امرناظ) تاتفصمل عللیدها ۲۵12017۵ 0ص 31و بر۳۷۷۵) آصماوهه ۳/2 11610 عقطا 0162060ض1 والناوع۲ مط 1 

۸۵۵۵۲۵ .2762 500160 مطا مر رعزظ) معصهاولوع؟ ممتاحامصمم متفصصتاوع م) وعاماهتنه۷ ۵11و ماماهاته27 ۲۵2011۷ 220 1۳00016206 

٩011 ۲ ۷۵۵6‏ ]0۲6010 10 ۵0615 ۸۵۵۲ ۵۶ برمه‌تتاهع2 مطا ,ولا ۵۶ ومبالم۷ ]عون مطا مه ۴۵۷۲9 ۵۶ ومباه۷ 10۷769۲ 6ظ) ها 
۰ :)جر وام0معه رن 4ص ]۱۷۲ صقط تعطعنط 


مولع 01۱ ,وماطامزنه۷ تدم رممنومی ماج تمطمن1]2 0۳2۷ رطمنعع1 100ظع40صع رصمتتمصصتا ۴۵ 1۵۲۱۷۵۲۵۵۰ 


نشریه آب و خاک 


ره 


(0 


هار۲ ساره که وهای سییر ۱۳۰۲ من ۲۶۹-۲۳ 


مدل‌سازی مقاومت فروروی خاک با استفاده از رگرسیون. شبکه عصبی 


۰ 


مصنوعی و برنامه‌ریزی بیان ژن 


شکراله اصفری ۱6*- مهسا حسنپور کاشانی (8- حسین شهاب آرخازلو 68" 


تاریخ دریافت: ۱۴۰۲/۱۱/۲۳ 


تاریخ پذیرش: ۱۴۰۲/۱۲/۲۹ 


چکیده 


اطلاع از مقاومت فروروی (۳1۶) خاک از نظر جوانه‌زنی بذر» رشد ريشه و عملیات خاک‌ورزی اهمیت فراوان دارد. اندازه‌گیری مستقیم ۳1 خاک 
به‌دلیل تغییرپذیری مکانی و زمانی شدید آن. کاری پرزحمت و گران می‌باشد. هدف از این پژوهش, ارایه مدل‌های رگرسیون خطی (۷115» شبکه عصبی 
مصنوعی (۸۱۷(۷) و برنامه‌ریزی بیان ژن (010۳)) برای برآورد ۳ از روی ویژگی‌های زودیافت خاک بود. در مجموع ۸۰ نمونه خاک سطحی (۰ تا 000 
۰) دست‌خورده و دستنخورده با استفاده از استواله‌های استیل به قطر و ارتفاع ۵ سانتی‌متر) از اراضی جنگلی, مرتعی و زراعی منطقه فندقلوی اردبیل در 
تابستان ۱۳۰۳ برداشته شد سپس برخی خصوصیات فیزیکی و شیمیایی زودیافت خاک در آنها اندازه‌گیری شد. مقاومت فروروی خاک به‌طور درجا در محل 
با استفاده از یک فروسنج مخروطی تعیین گردید و همزمان رطوبت خاک مزرعه در استوانه‌ها اندازه‌گیری شد. داده‌ها بهطور تصادفی به دو گروه آموزشی 
(۶۰ نمونه) و آزمونی (۲۰ نمونه) تقسیم گردید. مدل‌های 11۴ ۸۱۳ و 01۳ به‌ترتیب با استفاده از نرم‌افزارهای 5۳55 1۸۲۸۶ و عصعه 
195 26۳۲0 ایجاد شدند. نتایج مدل‌سازی نشان داد که رطوبت خاک مزرعه. سیلت و جرم مخصوص ظاهری نسبی» مهمترین متغیرهای ورودی در 
برآورد ۳1۳ خاک بودند. مقادیر آماره‌های ضریب تبیین (732), مجذور میانگین مربعات خطا (130۷1515» میانگین خطا (۵۷1۳) و ضریب نش-ساتکلیف (015) 
براساس داده‌های آزمونی برابر ۰/۴۴ ۷۲۳۵ ۸۱/۱۹ ۷۳۵ ۰/۱۹ و ۰/۳۶ ۰/۹۲ ۷۳۵[ ۰/۴۱ ۷۲۳۵ ۰/۰۵- و ۰/۹۲ ۰/۷۹ ۷۳۵ ۰/۹۱ ۱۷۳۵ ۰۱۱۳ 
و ۰/۶۳ به‌ترتیب برای بهترین مدل ]1۷۲ ۸۲ و طنآن) تعیین گردید. براساس نتایج آماره‌های ارزیابی مدل‌ها» می‌توان گفت که در منطقه مورد 
مطالعه» مدل ۸ از بیشترین دقت و مدل ,۷11 از کمترین دقت در برآورد +۳1 خاک برخوردار بود. 


واژه‌های کلیدی: تخمین. رطوبت جرمی, مقاومت خاک منطقه فندقلو, ویژگی زودیافت 


مقدمه 


برابر نفوذ یک فروسنج مخروطی یا مسطح نشان می‌دهد. از نظر جوانه 
زنی بذ رشد ريشه و عملیات خاک‌ورزی حائز اهمیت است. به‌طور 
کلی» اگر مقدار یک خاک از ۲/۵ مکاپاسکال. فراتر رود زشند و 


گسترش ریشه در خاک به‌طور معنی‌دار محدود خواهد شد ( ۱۷۷۵۱۱6 
7 ر .اه 61). 

مقاومت فروروی (۳۷۶) خاک با تغییر سایر ویژگی‌های خاک مانند 
کربن آلی. جرم مخصوص ظاهری و رطوبت خاک مزرعه» با زمان و 
مکان تغییر پیدا می‌کند. به‌طور مثال» اصغری و همکاران ( ۸52۱4 
17 ر ۵۶۵ گزارش گردند ضریب تفییرات (۷:)) مقاوست فروروش 


۱ و ۴ به‌تر تیب استاد و دانشیار گروه علوم و مهندسی خاک دانشکده کشاورزی و منابع طبیعی» دانشگاه محقق اردبیلی» اردبیل» ایران 


(:- نویسنده مسئول: 1۱۱۵6۴ ۵) زعقطع فقطه :اتفصط) 


۲- استادیار: گروه مهندسی آب دانشکده کشاورزی و منابع طبیعی دانشگاه محقق اردبیلی» اردبیل» ایران 
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4- ۳۵۵۵۲۲۵10۴ ۵۵۵ 


۳ شریه آب و خاک. جلد ۳۸. شماره ۲. خرداد - تیر ۱۴۰۳ 


در ۱۰۵ نمونه خاک برداشته‌شده (با فواصل ۲۰۰۰ متر) از اراضی زراعی 
دشت اردبیل (10 ۲۴ *< ۴۲) حدود ۵۰ درصد تعیین گردید و دامنه 
تغییرات کربن آلی از ۰/۳ تا ۲/۴ درصد جرم مخصوص ظاهری از ۱/۰۲ 
تا ۱/۶۳ گرم بر سانتی‌متر مکعب و ۳ از ۱/۱ تا ۶/۶ مگاپاسکال به 
دست آمد. همچنین اصغری و شهابی (2019 ,نطهطهط؟ ک نتهطعده) 
مقلدیر ۷ را برای مقاومت فروروی خاک در اراضی زراعی و بایر بخش 
شندآباد منطقه شبستر (10 ۴۰۰ < ۲۰۰۰) حدود ۴۶ درصد گزارش 
نمودند. از سوی دیگر اندازه‌گیری مستقیم ۳ نیز به‌دلیل خطاهای 
ابزاری» کاری طاقت‌فرسا و پرهزینه است؛ بنابراین استفاده از مدل‌های 
مختلف مانندرگرسیون خطی چندگانه (4(1/1.5 شبکه عصبی مصنوعی 
۸(۷((۳) و برنامه‌ریزی بیان ژن (012۳) برای برآورد 315 از طریق 
ویژگی‌های سهل‌الوصول و کم‌هزینه خاک می‌تواند در این زمینه 
راهگشا باشد. 

نتایج پژوهش بیات و همکاران (2008 ,۸ 6 2۷۵ظ) نشان داد 
که در برآورد ۳5 از روی تخلخل کل, درجه اشباع نسبی و رطوبت 
خاک مزرعه مدل ۸((۷ در مقایسه با مدل‌های رگرسیون خطی و 
غیرخطی از دقت بالاتری برخوردار بود. واز و همکاران ( .1 61 ۷22 
11 گزارش کردند که رطوبت خاک مهمترین ویژگی تأثیرگذار بر 
مقاومت فروروی خاک می‌باشد. سانتوس و همکاران ( .1 6 520105 
1 از جرم مخصوص ظاهری و رطوبت خاک مزرعه برای برآورد 
۴ خاک در مدل‌های رگرسیونی و ۸ استفاده کردند. آنها بیان 
کردند که ۳ با جرم مخصوص ظاهری رابطه مستقیم و با رطوبت 
خاک مزرعه رابطه معکوس دارد. شیخ‌زاده و همکاران ( ٩612206‏ 
6 .4 /6) گزارش کردند که در برآورد ۳1۶ خاک مقادیر ضریب 
تبیین (*8) و مجنور میانگین مربعات خطا (2151) به‌ترتیب در 
مدل ۸۱۲ برابر با ۰/٩۱‏ و 31۳۵ ۰/۳۷ و در مدل ۷1۴ برابر با ۰/۵۵ 
و ۷۳۵ ۰/۸۹ به‌دست آمد بنابراین مدل ۸۳ نسبت به مدل ۷11/۴ 
دارای دقت بالاتری بود؛ همچنین رطوبت خاک مزرعه» جرم مخصوص 
ظاهری و توزیع اندازه ذرات خاک مهمترین متغیرهای ورودی در برآورد 
۴ خاک تشخیص داده شد. اخیرا از سایر روش‌های هوشمند مانند 
الگوریتم ژنتیک و مدل 6012۳ در برآورد ویژگی‌های دیریافت خاک 
استفاده گردیده است؛ احمدی و همکاران (2016 ,اه ۶ ت4عصنطض) 
گزارش کردند که در برآورد هدایت هیدرولیکی اشباع خاک در دشت 
تبریزه مدل ۳) نسبت به مدل رگرسیون رایج به‌علت داشتن 2 بالا 
و 5۳ پایین دارای دقت بالاتری بود. آزاد مرد و همکاران 
(2019 ,.اه ۵ ۸220۳270) بیان کردند مدل ترکیبی ۸ با 
الگوریتم ژننیک در مقایسه با مدل ۷115 از دقت به مراتب بالاتری در 
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برآورد هدایت هیدرولیکی نزدیک اشباع خاک در دشت مغان استان 
اردبیل برخوردار بود. نتایج پژوهش احمدزاده کلیبر و فولاپناه 
(2023 مطفصدم‌ن۳12 عک تدماتهلمک1 طع02۵0مصطض) در دشت 
مقایسه با مدل‌های رگرسیون خطی و غیرخطی از دقت بالاتری در 
برآورد رطوبت‌های ظرفیت مزرعه و نقطه پژمردگی دائم برخوردار بود. 
زانگ و زانگ (2024 ,ع200 > ع22) گزارش کردند که مدل 
۳ توانست ضریب نفوذپذیری خاک را از روی برخی ویژگی‌های 
خاک مانند نست پوکی» اندازه موّثر و میانگین اندازه ذره با دقت بالاتری 
9 برابر با ۰۴« 9 ۳۰(۷9۲ پرابر با ۰۰۴۸ 0/۰( نسبت به سایر مدل‌های 
هوشمند و روابط ریاضی تخمین بزند. 

آگرچه مدل‌های رگرسیونی و شبکه عصبی مصنوعی (۸(۷۵) 
مختلفی برای برآورد ۳۶ خاک در نقاط مختلف دنیا ارائه گردیده است 
ولی کاربرد عملی این مدل‌ها در خاک‌های سایرمناطق جهان با شرایط 
اقلیمی متنوع» همواره توسط محققان با شک و تردید همراه بوده است. 
از سوی دیگر تحقیقات انجام گرفته در زمینه استفاده از مدل برنامه 
ریزی بیان رن (0۳۲۳۲) برای برآورد ۶ خاک انگشت‌شمار است. 
بررسی منابع نشان داد که تاکنون تحقیقی در زمینه برآورد ۳5 توسط 
مدل‌های هوشمند مانند ۸ و ۳۲۳ و مقایسه آن با مدل‌های 
رگرسیونی در جنگل‌های فندقلوی اردبیل اجرا نگردیده است لذا موضوع 
پژوهش حاضر از این نظر دارای نوآوری است. با توجه به اينکه منطقه 
فندقلوی استان اردبیل از مناطق گردشگری مهم و پرتردد کشور بهویژه 
در فصل تابستان محسوب می‌شود و تغییر کاربری اراضی جنگلی و 
مرتعی به زراعی نیز در آن با شدت زیاد در حال وقوع است لذا ارائه 
مدل‌های هوشمند برای اطلاع از وضعیت فیزیکی و مکانیکی خاک‌های 
منطقه مذکور حائز اهمیت فراوان است. بنابراین اهداف پژوهش حاضر 
عبارت بودند از (۱) ارائه مدل‌های 0۷11/۶ ۵۷۲ و طت) برای برآورد 
۴ از روی متفیرهای زودیافت در خاک‌های جنگلی» مرتعی و زراعی 
منطقه فندقلوی استان اردبیل (۲) مقایسه دقت مدل‌های مذکور در 
برآورد ۳۶ خاک با استفاده از معیارهای 127 ۷5۳( میانگین خطا 


مواد و روش‌ها 
در شرق استان اردبیل (مختصات جغرافیایی ۴۵۳ ۲۲ ۴۸ تا ۵۳ ۲۳ 
۴ طول شرقی و ۱۰۳ ۲۴ ۲۸۶ تا "۲۵ ۲۴ ۲۸۳ عرض شمالی) قرار 
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اصغری و همکاران» مدل‌سازی مقاومت فروروی خاک با استفاده از رگرسیون. شبکه عصبی مصنوعی و برنامه‌ریزی بیان ژن ۲۷۳ 


دارد که فاصله آن از مرکز استان اردبیل در حدود ۳۰ کیلومتر می‌باشد. 
تعداد ۸۰ نمونه خاک دست‌خورده (با حداقل دست‌خوردگی) و دست 
نخورده با استفاده از استوانه‌های استیل به قطر و ارتفاع ۵ سانتی‌متر) 
از محدوده عمق ۰ تا ۱۰ سانتی‌متر اراضی شیبدار و تپه‌ای با کاربری 
جنگلی (۲۰ نمونه)» مرتمی (۳۳ نمونه) و زراعی (۳۷ نمونه) به هم 
چسپیده به فواصل حدود ۵۰ متر و به طور تصادفی- منظم در تابستان 
۲ برداشته‌شد. موقعیت هر یک از نقاط نمونه‌برداری‌شده با استفاده 
از دستگاه موقعیت‌یاب جهانی (0۳5)) ثبت گردید (شکل ۱). 

پس از هواخشک کردن نمونه‌های خاک بخشی از آن‌ها برای 
تعیین پایداری خاکدانه‌ها جدا گردیده و بخش دیگر برای تعیین سایر 
ویژگی‌های فیزیکی و شیمیایی از الک 91 ۲ گذرانده شد. توزیع اندازه 
ذرات شن, سیلت و رس و کلاس بافت خاک به روش هیدرومتری۴ 
قرائته (2002 ,0 6 66)) تعیین و پارامترهای میانگین هندسی قطر 
(یل) و انحراف استاندارد هندسی (یه) ذرات خاک از روی درصدهای 
شن» سیلت و رس محاسبه گردید (1985 ,العاحصصهن)). کرین ۳ به 
روش اکسیداسیون تر (1982 ٩010۳0675,‏ عک ۵1500)» اسیدیته (0۳۲) 
و هدایت الکتریکی (۳0) در عصاره ۱ به ۲/۵ خاک به آب و آهک به 
روش تیتراسیون با اسید و سود (1985 ,۳۵۵6) تعیین شد. جرم 
مخصوص ظاهری («1ظ) در استوانه‌های استیل ( ,۲۲۵۲/86 6 1216ظ 
۳۶ و جرم مخصوص حقیقی (۳10) به روش پیکنومتر ( > 1216 
19960 ,26 ت۳2) اندازه‌گیری شد. جرم مخصوص ظاهری نسبی 
([12:6ظ) از روی داده‌های جرم مخصوص ظاهری و درصد رس 
محاسبه گردید (2014 ,. ۶ ن24حاحتاعک). میانگین وزنی قطر 
(1۷]۷۷۲) و میانگین هندسی قطر (01۷۲1) خاکدانه‌ها به روش الک 
خشک در ۵۰ گرم خاک هواخشک گذرانده شده از الک ۴/۷۵ میلی‌متر 
با استفاده از دستگاه شیکر الک با به کارگیری الک‌های با قطر سوراخ 
و محاسبه گردید (1936 ,۷۵067). مقاومت فروروی (۳3) خاک به‌طور 
مستقیم در مزرعه با استفاده از دستگاه فروسنج مخروطی مدل عقربه‌ای 
ساخت شرکت اکل کمپ کشور هلند در ۵ تکرار اندازه‌گیری شد 
(2002 ,طمو:۱۷]0 ی تع/«م۲]) و همزمان رطوبت خاک مزرعه در 
حین اندازه‌گیری با استفاده از استوانه‌های خاک برداشته شده برای 
اندازه‌گیری 31 پس از آون خشک کردن در دمای 0 ۵ به روش 
وزنی تعیین گردید (1986 ,01270067)). 

نرمال بودن توزیع داده‌ها با آزمون کولموگروف- اسمیرنوف بررسی 
گردید سپس محاسبات مربوط به آمار توصیفی و ضریب همبستگی 
پیرسون بین متغیرها و نیز اشتقاق مدل‌های رگرسیونی خطی چندگانه 
(1]5) بین متغیرهای زودیافت خاک با متغیر دیریافت (۳1۴) به روش 
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گام به گام در نرم‌افزار 22 5۳55 انجام شد. 

مدل‌های شبکه عصبی مصنوعی (۸(0 با استفاده از نم‌افزار 
۶ ایجاد گردید. در این پژوهش شبکه پرسپترون سه لایه 
ام یک لاه ورودی یک لاه هام و یک لاه شروش زا لگوریتج 
پس انتشار خطا و ۱۰۰۰ تکرار برای رسیدن به حداقل خطا و حداکثر 
پایداری شبکه به کار گرفته‌شد و تعداد نرون‌های مناسب به روش سعی 
و خطا در لایه پنهان برای سه مدل شبکه عصبی به‌ترتیب برابر ۶ ۵ و 
۲ انتخاب گردید. نرون‌های پنهان داخل شبکه پرسپترون چند لایه, در 
هر لایه به نرون‌های لایه بعد متصل هستند. به‌منظور تجزیه و تحلیل 
ویژگی‌ها در داده‌های ورودی» توابع فعال‌سازی تانژانت سیگموئید و 
خطی به‌ترتیب در لایه‌های پنهان و خروجی هر سه مدل به‌کار برده شد 
(2019 ,.۵1 ۵۶ نصدطعمطن). به‌منظور آموزش شبکه با دقت بالا از 
الگوریتم آموزشی مار کوآرت-لونبرگ استفاده گردید. 

برنامه‌ریزی بیان ژن (0۳۳)) توسط فریرا (2001 بدتذع۳۵) 
پیشنهاد شد. این روش یک الگوریتم توسعهیافته براساس برنامه‌ریزی 
ژنتیک و الگوریتم ژنتیک می‌باشد و مبنای تمامی آنها براساس نظریه 
تکامل داروین استوار است که در آن کروموزوم‌های خطی و ساده با 
طول ابت مشابه آنجه در الگوریتم ژننیک استفاده می‌شود و ساختارهای 
شاخه‌ای با اندازه‌ها و اشکال متفاوت مشابه با برنامه‌ریزی ژنتیک ترکیب 
می‌شوند. برای مطالعه بیشتر در زمینه جزئیات مدل ۳۳) می‌توان به 
مقاله زانگ و زانگ (2024 ,ع22 > 2۳۵۵) مراجعه نمود. در این 
تحقیق برای مدل‌سازی مقاومت فروروی (۳۷۶) خاک براساس برنامه 
ریزی بیان ژن (0۳۳) از نرم‌افزار 4.0 10015 2070 0626 استفاده 
شد. برای این منظور از مجموعه عملگرهای ریاضی شامل چهار عمل 
اصلی (+۰ - * و /» جذر اکسپونشیال لگاریتم طبیعی» توان ۲ و ۲ و 
ريشه سوم استفاده شد. همچنین پارامترهای لازم برای مدل‌سازی 
۳ در جدول ۱ آورده شده است. 

در همه مدل‌های به‌کار رفته در تحقیق حاضر از مجموع ۸۰ داده 
رای هر متفیر, بعد از تصادفی کردن کل داد‌ها در نم‌افزار اکسل؛ ۶۰ 
داده برای آموزش مدل و ۲۰ داده برای آزمون مدل در نظر گرفته شد. 
برای ارزیابی دقت و کارآیی مدل‌های ۸۷11 ۸۸۷۷ و 017۳ در برآورد 
خطا (80050), مينکین خطا (۸415) و ضریب نش-ساتکلیف (۱05 
(2019 .41 6۲ نصدح:0)) استفاده شد. مقادیر بالای 187 (نزدیک به 
۱( بیانگر دقت بالای برآورد می‌باشد. 
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شکل ۱- موقعیت جغرافیایی و نقشه توزیع نقاط نمونه‌برداری در تصویر ماهواره‌ای منطقه مورد مطالعه 
2 5000160 6ظ) ۵۶ 6عمصصز ماتام)هو صا امن 60وی ۵ مرح ممتابطاتتاعنل 0صه صمتاتعمن منصامرمنعه۵ع) -1 ۳۲6اع۲1 


جدول ۱- مقادیر پارامترهای مدل برنامه‌ریزی بیان زن (01۳) 
م۵2۱۵ ۳۵061 ای ۵ ومبتاه۷ 106 -1 12۳016 


تنظیمات عمومی مقادیر 

اصمصصامن ز20 موی ۷۵6( 
تعداد کروموزم‌ها 90۴6محصمتط ۵ تعطاصصت لا 30 
اندازه سر 5126 ۲1620 8 
تعداد ژن‌ها 6065ع 0۶ تعطاصصن( 3 

تابع پیوند طمتاعصا عصل01] جمع (+) 
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از لحاظ آماری کمترین مقدار 3۷۲8۳ صفر می‌باشد که در آن 
وضعیت مقادیر اندازه‌گیری و برآورد شده ۳18 باهم برابر هستند؛ هرچه 
مقدار 1۷157 کوچکتر باشد برآوردهای انجام شده توسط مدل برآورد 
کننده دقیق‌تر است. مقادیر مثبت 1۷1۳ بیش‌برآوردی و مقادیر منفی آن 
کم‌برآوردی مدل را نشان می‌دهد. اگر ضریب 5 برابر یک باشد بیانگر 
آن است که مدل‌سازی با دقت بالایی انجام گرفته و مقادیر اندازه‌گیری 
و برآورد شده ۳ در تطابق کامل هستند. مقدار برابر با صفر 5 
نشان‌دهنده آن است که برآوردهای مدل» دقتی برابر با دقت ميانگین 
داده‌های اندازه‌گیری شده دارند. 
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عوامل ژنتیک مقادیر 
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نرخ جهش ۲۵16 1۷۲۵۵1100 00094 
نرخ وارون‌سازی ۲۵/6 12۷65100 01 
نرخ ترکیب تک نقطه‌ای ۲۵06 ۲660002/00 اصزهم عم 03 
نرخ ترکیب دو نقطه‌ای ۲۵06 طمنامهزطعع۲ ]ص00 11۷/0 03 
نرخ ترکیب ژن ۲۵/6 حمنادصنطهعع۲ عع0 01 

1 
هن 7 (زلا ۳3 
)۲( ۱۷۲۳ 
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رس ۷ - رس ۷) و 5 
رس رل 
در معادله‌های ۱ ۲ و ۲ :۰۷ ری:۷ و لا به ترتیب مقادیر اندازه 
گیری» برآورد شده و میانگین مقادیر اندازه گیری شده ۳ توسط مدل 
و 1 تعداد داده‌ها برای ۶ در هر مدل می‌باشد. از میان مدل‌های ایجاد 
شده برای برآورد ۳ مناسب‌ترین مدل 211۴ ۵۱۲ و طتلن با 


1 < ول 


۳) 


ترکیب متفیرهای ورودی در داده‌های آزمونی براساس محک‌های 
آماری مذکور تعیین گردید. آزمون مدل‌ها در تشخیص دقت مدل به 
ازای داده‌هایی که در آموزش مدل استفاده نشده‌اند به کار می‌رود. 
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آماره‌های توصیفی ویژگی‌های خاک 

توصیف آماری ویژگی‌های خاک منطقه فندقلوی اردبیل در جدول 
۲ نشان می‌دهد که داده‌های آموزشی و آزمونی به کار رفته در توسعه و 
ارزیابی مدل‌های رگرسیون خطی چندگانه (۷1۶» شبکه عصبی 
مصنوعی (۸۷۱۷) و برامه‌ریزی بیان ژن (02۴) درای توزیع نسبتاً 
مناسبی هستند و این موضوع از نظر توسعه و ارزیابی مدل‌ها حائز 
اهمیت است؛ به‌طور مثال. دامنه تغییرات بالا به‌ترتیب در داده‌های 
آموزشی و آزمونی برای مقاومت فروروی خاک (۳) از ۱/۲ تا ۴/۸ 
کال از ۳ ۵ عکارایکال رارقا رتخا 
(۳۷0) از ۱۱/۴۸ تا ۵۰/۴۴ درصد وزنی و ۱۲/۰ تا ۵۰/۴۷ درصد 
وزنی و کربن آلی خاک از ۱/۰۱ تا ۷/۱۷ و ۱/۱۵ تا ۶/۲۱ درصد موّید 
این موضوع است. همچنین مقادیر ميانگین داده‌ها در سری آموزش و 
آزمون به هم نزدیک بوده و به‌ترتیب برای ۳1 ۲/۱۵ و ۲/۲۲ 
مگاپاسکال, برای 1۳ ۲۳/۴۳۶ و ۲۵/۷۳ درصد وزنی و برای کرین 
آلی. ۲/۸۵ و ۲/۹۰ درصد می‌باشد. لازم به توضیح است که مقادیر 
بالای کربن آلی به کاربری جنگلی و مقادیر پایین کربن آلی به کاربری 
زراعی مربوط می‌باشد. همچنین مقادیر بالای به‌دست آمده برای 


رطوبت خاک مزرعه (در لحظه اندزه‌گیری ۳) ب‌علت وقوع بارندگی 
در اواخر تیر ۱۳۰۲ در منطقه مورد مطالعه می‌باشد. مقادیر متوسط 
میانگین وزنی (۱81۷19) و هندسی (011) قطر خاکدانههابه‌ترتیب 
در داده‌های سری آموزش برابر با ۱/۵۵ و ححط ۱/۰۳ و در داده‌های 
سری آزمون برابر با ۱/۶۱ و حححط ۱/۴ به‌دست آمد (جدول ۲)؛ مقادیر 
بالای شاخص‌های پایداری خاکدانه‌ها (۷۲۷۷۲ و 01۷۲۲) در تحقیق 
حاضر در مقایسه با سایر خاک‌های استان اردبیل را می‌توان به کربن 
آلی زیاد خاک در منطقه مورد مطالعه نسبت داد. شیخ‌زاده و همکاران 
(2016 .۵1 ۶ 12۵06ع) میانگین کربن آلی و ۳ خاک را 
به‌ترتیب ۱/۱۱ درصد و ۲/۷۸ مگایاسگال در ۱۰۵ نمونه خاک زراعی 
برداشته‌شده از دشت اردبیل گزارش کردند. امیرعابدی و همکاران 
(2016 .۵1 ۶۶ ن0عطهتن۸) نیز ميانگین کربن آلی و ۷۱۷ را به 
ترتیب ۰/۷۷ درصد و 10۳0 ۰/۵۰ در ۱۰۰ نمونه خاک زراعی برداشته 
شده از دشت اردبیل گزارش کردند. 

با توجه به شکل ۲ مشاهده می‌شود که کلاس بافت خاک‌های 
مورد مطالعه براساس مثلت بافت آمریکایی لوم» لوم رسی سیلتی لوم 


شنی لوم رسی» لوم سیلتی و رسی سیلتی بود که بیانگر توزیع گسترده 
بافتی در منطقه مورد مطالعه است. 


جدول ۲- آماره‌های توصیفی متغیرهای خاک در منطقه مورد مطالعه 
مسج 0عنباو عمط رز معاطاهتیه۱ م5 ۵۶ متاوناهاه 00۷۵ز۵عع(1 -2 12016 


آزمون 


عصنامع ]1 

90۰ ارگ ۱/۸ ۱/۳۹۹ ۱/۱ 0۰ 
و12 25.73 320.47 12.90 1024 
010 063 079 043 0017 
017 241 208 2.02 015 
1127 2028 38.40 2394 9353 
1235 45 146 25-3 12.90 
1189 5/۸ 6480 15-31 13-52 
005 006 021 002 005 
3220 1175 2043 6094 3229 
025 11 198 119 028 
01 104 12 08 011 
171 290 631 115 175 
2017 614 15.6 21 116 
043 633 1.34 58 038 
00 011 033 004 019 
152 222 53 1.03 085 


ی ِِ ۲۷ 
اااگ ۱/۸ ۱/۳۹ ۱/۱ تسه 
2146 5044 148 (1ع ع ۳۷۷۵۲۹ 

063 015 044 9/۵ 
244 201 201 (تصه ع) ۳۲ 
2118 2-18( 295 (۶0) «012) 
)٩6( 21-1 6519 4462‏ ]11 
2339 6410 1314 (۶6) 5۵00 
006 018 001 (صص) ی 
1175 203 694 () یه 
135 2:31 09 (صصص) ۷۲۷۷۲۵ 
103 133 08 (صصه) ن(0 
285 1-17 1.01 050 0 
1.18 1032 32-12 (60) :0200 
629 12 545 ۱ 
015 108 002 (ه 45) ومن۴ 
215 18 12 (۷۲۳۵) ۳ 


61 رطوبت خاک مزرعه (احعاطه ۷/2۲ 1610), 1ع:۳1:جرم مخصوص ظاهری نسبی (0605107 عل ۲61201۷6م(]: جرم مخصوص حقیقی 


(راتعجعل ع1علایدع)» ملن: رس ]911: سیلت» ۹200: شن, ول: میانگین هندسی قطر ذرات (7عاعه1 حدعحط 16تع0۳0عع)» وه: انحراف استاندارد هندسی 
ذرات (06۷1210108 4تع‌صهاه عتتاع0عع)» ۷۲۷۷۲۵ میانگین وزنی قطر خاکدانه‌ها به روش الک خشک (012۳06167 )ظع۷61 صدع۳)» ]0۷1: ميانگین 
هندسی قطر خاکدانه‌ها (ععاعصونل صععط عنتاعصصمعع). :06: کرین آلی (طهطنعء عنصعععم و20/0): کربنات کلسیم معادل» *۳1: مقاومت فروروی 
(ععصماعلدع۲ طمتاح0اه0ع0). 90: انحراف استاندارد. ۷620 میانگین, ۷2 حداکتر ۷۲18 حداقل. 
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دما 0صج٩‏ 
‌‌ ۰ 
لوم سئی 


5290 )۶۰( 


شکل ۲- توزیع کلاس بافت در خاک‌های مورد مطالعه 
۶06 500160 66 جر دعقل آهساه) ۵۴ جمتاتطاتت)6 -2 ۲12۱۲۵ 


ضرایب همبستگی پیرسون (7) بین متفیرهای خاک در جدول ۳ 
نشان می‌دهد که همبستگی خطی قوی (۵.60< ۲) بین شن با سیلت 
و میانگین هندسی قطر ذرات (ی) سیلت با انحراف استاندارد هندسی 
متغیرها به‌عنوان ورودی مدل به‌طور همزمان برای ایجاد مدل‌های 
ضرایب برآوردی تا چه حد نسبت به‌حالتی‌که متفیرهای برآوردی 
همبستگی خطی ندارند افزايش یافته است (2014 ,.اه 6 ۷2270121 ). 
نتایج آنالیز هم‌خطی حداکثر تورم واریانس (۷1۳) مدل‌های 1/۲ 
مربوطه (جدول ۴) مقادیر کمتر از ۵ را نشان دادند که مقدار حداکثر آن 
در این پژوهش ۱/۴ به‌دست آمد. به‌طور کلی اگر آماره ۷1۳ برای یک 
متغیر کمتر از ۵ باشد آن متغیر در مدل وارد می‌گردد در غیر این صورت 
باید از مدل حذف شود. 


۱۳۹ 


نتایج مدل‌های رگرسیون خطی چندگانه (۷118) برای 
برآورد مقاومت فروروی (۳1۴) خاک 

به‌طور کلی براساس روش گام به گام در رگرسیون خطی» ۳ نوع 
مدل ۷ برای برآورد ۴ با استفاده از متغیرهای رطوبت خاک 
مزرعه (4۳0#6 سیلت و جرم مخصوص اهری نسبی (8۳201) از 
بین متغیرهای زودیافت در منطقه مورد مطالعه تشخیص داده شد 
(جدول ۴). نتايج ضریب همبستگی پیرسون بین ۳ با متغیرهای 
و معنی‌دار بین ۳۶ با ۳۷۷ (**۰/۴۳۵- - ) یافت شد. همچنین 
همبستگی منفی و معنی‌دار بین ۳ با سیلت (**۰/۲۶- < ۲ و 
همبستگی مثبت و معنی‌دار بین ۳۶ با انحراف استاندارد هندسی ذرات 
(وه) (**۰/۳۶ 2 ) به‌دست آمد؛ البته به‌دلیل همراستایی سیلت و 6 
از یه به‌عنوان ورودی در تخمین ۳ استفاده نشد و سایر متغیرهای 


(۷]۳) 120001 طمتاحاگص1 معصهاجه ۷ -2 
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ورودی (به استقنای آه»(181) نیز براساس تشخیص روش کام به کام و 
احتمالا به علت نداشتن همبستگی معنی‌دار با ۳1۶ در مدل‌های ]1۷۲ 
ظاهر نشدند. 

مشاهده می‌شود ۳۷۷ در هر سه مدل ]۷۲ به‌عنوان متغیر 
ورودی و با علامت منفی ظاهر گردیده است (جدول ۴) بنابراین در 
پژوهش حاضر رطوبت خاک مزرعه (۳1۷۲) مهمترین عامل اثرگذار در 
برآورد مقاومت فروروی خاک بوده است. در مدل‌های رگرسیونی ارائه 
شده توسط بیات و همکاران (2008 ,. 6 2۵۲ واز و همکاران 
(2011 .اه 6 ۷2) و شیخ‌زاده و همکاران ( .1 ۶ 12۵060«عط٩‏ 
06 نیز رطوبت خاک مزرعه یکی متغیرهای ورودی مهم در برآورد 
۴ گزارش گردیده است. علامت منفی ۳۷۷ در هر سه مدل .]1۷۲ 
نیز بیانگر آن است که با افزایش رطوبت خاک محل و نرم شدن خاک» 
مقاومت مکانیکی خاک در برابر نفوذ فروسنج مخروطی کاهش یافته 
است که با یافته‌های سانتوس و همکاران (2011 ,.۵1 21 920105) و 
اصغری و شهابی (2019 ,1طحطدط٩‏ > تتدطعمض) در یک راستا می 
باشد. همچنین مشاهده می‌شود در مدل .1۷11 شماره ۲ و ۲ سیلت با 
علامت منفی به‌عنوان متغیر ورودی وارد مدل گردیده است (جدول ۴)؛ 
ذرات سیلت در مقایسه با زبری ذرات شن و چسپندگی ذرات رس از 
خاصیت نرمی برخوردار هستند (1985 ,02۳0061)) بنابراین انتظار 


۳۷۳۷ 


می‌رود با افزایش درصد سیلت در خاک به علت ایجاد حالت سست و 
نرم مقاومت فروروی خاک کاهش يابد. در مدل .1۷11 شماره ۲ علاوه 
بر :۷6 و سیلت» جرم مخصوص ظاهری نسبی (۲31«:61) نیز با 
علامت منفی وارد مدل گردید. به‌نظر می‌رسد با افزایش سهم رس در 
توزیع اندازه ذرات خاک جرم مخصوص ظاهری مرجع در مخرج کسر 
محاسبه 31:61 کاهش و در نتیجه 317701 افزایش و منحر به کاهش 
۶ می گردد. کلیشادی و همکاران (2014 ,.41 6 1150201ع1) توضیح 
نسبتاً کاملی درخصوص کاستی‌های جرم مخصوص ظاهری و مزیت 
های استفاده از 11:01 به‌جای آن ارائه کرده‌اند. همان‌طوری‌که از 
جداول ۳ و ۴ نمایان است برخلاف انتظار, در این پژوهش بین کربن 
آلی و ۳ همبستگی معنی‌دار یافت نشد بنابراین کربن آلی نتوانست 
به‌عنوان یک متغیر ورودی مدل‌های رگرسیونی تخمین‌زننده ۳15 ظاهر 
گردد که دلیل اين مسأله را شاید بتوان به ضریب تغییرات بالای آن 
(6070 <۷)) در سه کاربری مجاور هم جنگلی» مرتعی و زارعی منطقه 
فندقلوی اردبیل نسبت داد. این در حالی است که در تحقیق شیخزاده و 
همکاران (2016 ,.لة 66 9۳6۳2۵06۳) در خاک‌های دشت اردبیل» 
بین کربن آلی و 5[ همبستگی منفی و معنی‌دار (*۰/۴۷۴- <۲) پیدا 
شد و کربن آلی در برخی مدل‌های ۱۷1 برای برآورد ۳15 به‌عنوان 
ورودی مدل انتخاب گردید. 


جدول ۳- ضریب همبستگی پیرسون )۳ بین متغیرهای مورد مطالعه (۵-80) 
(۵<80) ۲۵۵0۱6 0۱60 م۵۲ رح اصعتتلهی صمتاهاه‌نی ممعروع۳ -3 12016 


1۳۳۷0 (۵۱ ۳1 [۷۲۷۷۲ 0] و2060‎ 
-0.45  . 7 019 00 0.01 -3 
0078۳ 050 05 5 014 040 
0.44" .9 0.27 018 011 042 
-7 -7 0.07 0.24 -0.27* ۰. 8 
0597 015 004 018 -029* 
0.47 0۵4 0۵6 2 -3 045" 
-012 033" 0.05 034 0.36" ۰.5 

1 0.577 37 0 017 0037" 
1 0.527 3 001 -9 
1 014 008 -.8 
1 0.877 
1 002 
1 


# 4 191 11 ا ۱ 
036 010 014 6 ۰-1 0.06- ۳۴ 
3 "0.55 *9- 4 0.597 1 0 
2 0.887 "6 *6.- 1 9200 
4 0.50 1-- 1 11 
"*(- 1 0۳7 
001 1 1 
1 6 
1۳۷۷ 
متا 
(۳1 
1-۳۳۳۵ 
00 
0200 


#و ##: به‌ترتیب معنی دار در سطح احتمال ۱ و ۵ درصد. توصیف علائم متغیرها در زیر جدول ۲ آورده شده است. 


۰ ,0.01 > ظ 20 0.05 > ظ 24 اصهمنصعز٩‏ : 


۰ .2 12016 ۲مصنا ماطهاته2۳ ود صمتام‌تهوع عامطصیره عاصاهتته ۱۷ 


جدول - مدل‌های رگرسیون خطی چندگانه (۲1.18) مقاومت فروروی (۳18) خاک 
۴ 01 0۴ ۴0۵0615 ,]۱۷۲ -4 12016 


حداکثر ۷۲۳ مدل شماره 
۲ *۷]2۵ ۷۲001 مداصصت ۱ 
1.00 0 0.0355 - 3.007 ۳۳2 1 
103 از 0.024 ۳۷۷ 0.041 - 4.22 ۳2 2 
140 -8.161 -۳۳ ۲12761 5.098 - 511 0.030 - ۳۷۷/۵ 0.059 3 


توصیف متغیرها در زیر جدول ۲ آورده شده است. 


۰ 12016 ۱۱06۲ م2721 وا ممتام‌تیعععل ععاماهتته ۲۷ 
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جدول ۵- معیارهای ارزیابی مدل‌های ر گرسیونی خطی چندگانه (۷11.18) مقاومت فروروی (۳1) خاک 
۴ 6 ۴0۵061 ۱۷۲۲/۱۸ ۲0۲ هنن صمتامتاه۳۲۵ -5 12016 


آموزش آزمون 
3 
۳ میم هم ۲ کا ‏ ومیم 
018 0001 077 8 020 .۰۰ 0.12 
031 0003 030 31 036 ۰ 019 
045 0 03 5 029 019 


1690102 6 


۹ تچ متغیرهای ورودی شماره 
(0۳۴۵ معاطاهننه۱ )رح صاصصت ۱ 
132 027 0/۵ 1 
119 044 1 ,)۳۷۷ 2 
125 034 اعساظ باه ,۲۱۷۷۵ 3 


5 ضریب نش-ساتکلیف, 1[۷]5۳: مجذور ميانگین مربعات خطاء ۷1۳ میانگین خطاء 7: ضریب تبیین. توصیف متغیرها در زیر جدول ۲ آورده شده است. 
معاطافتته ۷ .اممتمتگگهمه ککنامانک مه :۱۱5 ختمعه صفمصه :ظا زرمعه متمناوو صفعصط اممر ,5۳۲5۴ زصمتاحصتصماعل گ اصمئ‌تگلهمع :۳8۰ 
۰ ۲2016 06۲صنا ماطاهاته2۷ ود متامت0وع0 


در جدول ۵ مقادیر معیارهای ضریب تبیین (4157 مجذور میانگین 
مربعات خطا (18151) میانگین خطا (۷1۳) و ضریب نش- ساتکلیف 
(5() برای ارزیابی مدل‌های ۷11 از نظر دقت برآورد 18 آورده شده 
است. استفاده از 182 به تنهایی برای سنجش اعتبار مدل‌ها ممکن است 
فریبنده باشد و الزاماً 182 نزدیک به ۱ دلیل بر تطابتق کامل مقادیر 
اندازه‌گیری با مقادیر برآوردی ۳1 نیست یعنی بالا بودن 13 دلیلی بر 
کم بودن خطا نیست (1996 .617 1>021) بنابراین در پژوهش حاضر 
رتبه‌بندی مدل‌های *211 ۸۵۸۸ و 0 در برآورد ۳8 بر مبنای 
پایین بودن 1۷15۳ و بالا بودن ۷5 انجام شد. براساس پایین بودن 
میزان 1۷15 و بالا بودن میزان 5 در داده‌های آزمونی» مشاهده 
می‌شود مدل ۱۸18 شماره ۲ (جدول ۵) با متفیرهای ورودی رطوبت 
خاک مزرعه و سیلت از نظر برآورد ۳18 دارای دقت بالاتری در مقایسه 
با سایر مدل‌های .۷11 بود. مقادیر 122 1216 21۳ و 5 این مدل 
برای داده‌های آموزشی به‌ترتیب برابر ۰/۳۱ ۷۳۵ ۰/۷۰ ۷۲۳۵[ ۰/۰۰۳ 
و ۰/۳۱ و برای داده‌های آزمونی برابر ۰/۴۴ ۷۳۵ ۸/۱۹ 1۷۲۳۵ ۰/۱۹ 
و ۰/۲۶ تعیین گردید (جدول ۵). مقادیر مثبت 1۷1۳ برای داده‌های 
آموزشی و آزمونی بیان‌گر برآورد بیشتر از مقدار واقعی ۳1 خاک توسط 
مدل شماره ۲ می‌باشد. همچنین معادله رگرسیونی این مدل (جدول ۳) 
نشان می‌دهد که تأثیر منفی متفیر رطوبت خاک مزرعه به علت داشتن 
ضریب رگرسیونی بالاتر در برآورد 315 به مراتب بیشتر از متغیر سیلت 
می‌باشد. در رتبه‌های بعدی» مدل شماره ۳ با متغیرهای ورودی رطوبت 
خاک مزرعه, سیلت و جرم مخصوص ظاهری نسبی و مدل ۱ با متفیر 
ورودی رطوبت خاک مزرعه قرار دارد (جدول ۵؛ پس در مجموع 
براساس آماره‌های ارزیابی سه مدل رگرسیونی فوق‌الذکر مدل شماره 
۲ برای برآورد ۳۷۶ ارجح‌تر است. در پژوهش شیخ‌زاده و همکاران 
(2016 ,۵1 ۰۶ طع0همطو5۳0) نیز بهترین مدل ۱118 با متغیرهای 
ورودی رطوبت خاک مزرعه و تخلخل کل با مقادیر 57 1۷۲5۲ و 
5 به‌ترتیب برابر با ۰۰/۴۱ ۷۳۵ ۰۱/۱۶ 1۳۵ ۰/۱۷ (داده‌های 


آزمونی) برای برآورد ۳18 در خاک‌های زراعی دشت اردبیل تعیین گردید. 


نتایج مدل‌های شبکه عصبی مصنوعی (۸(۷۲) برای برآورد 
مقاومت (۳10) خاک 

برای ایجاد مدل‌های شبکه عصبی مصنوعی (۲(آ(ظ) برآورد کننده 
از ۳ مدل با همان متغیرهای ورودی مدل‌های ,۱۷11 (جدول ۵) 
استفاده شد؛ جدول ۶ معیارهای ارزیابی ۲ مدل ۸۸۲ ارائه شده برای 
برآورد ۳ را نشان می‌دهد. رتبه‌بندی مدل‌های ۸۷ نیز همانند 
مدل‌های 11.1 (جدول ۵ بر مبنای پایین بودن 1157 و بالا بودن 
5 انجام شد. براساس معیارهای فوق در داده‌های آزمونی (جدول ع۶)» 
بهترین مدل ((۸ برای برآورد ۳18 مدل ۸۱۱ شماره ۳ با متفیرهای 
ورودی رطوبت خاک مزرعه, سیلت و جرم مخصوص ظاهری نسبی در 
مقایسه با سایر مدل‌های ۸۱۱ بود. مقادیر 182 121517 ۷15 و 5 
این مدل برای داده‌های آموزشی به‌ترتیب برابر ۰/۷۹ ۷۳۵[ ۰۱۲۹ 
۵ ۰/۰۵ و ۰/۷۸ و برای داده‌های آزمونی برابر ۰/۹۲ ۷۳۵ ۰۱۴۱ 
52 ۰/۰۵ - و ۰/۹۲ تعیین گردید (جدول ۶). مقادیر مثبت و منفی 
۲ برای داده‌های آموزشی و آزمونی به‌ترتیب بیان گر برآورد جزتی 
پیشفر و کت از مقتاز واقتی *ظ خاک بط من شماره ۴ م‌بانی: 
در رتبه‌های بعدی به‌ترتیب مدل ۸٩‏ شماره ۲ و شماره ۱ دارای 
دقت بیشتری در برآورد *۳3 می‌باشد؛ مدل ۸۲ شماره ۲ نشان می 
دهد اضافه شدن متغیر سیلت به مدل ۸۸۲ شماره ۱ موجب کاهش 
خطای مدل و افزایش دقت آن گردید به‌طوری که براساس داده‌های 
آزمونی» 1151۳ مدل (۸(۱ شماره ۲ نسبت به مدل (۸(۱ شماره ۱ 
به میزان ۴۵ درصد کاهش و ضریب "٩‏ مدل ۸۲ شماره ۲ نسبت 
به مدل ۸ شماره ۱ به میزان ٩۳‏ درصد افزایش یافت. بیات و 
همکاران (2008 .21 6 2۵۲ظ) گزارش کردند افزودن متغیر درجه 
اشباع نسبی به مدل ۸۱ برای برآورد ۳3۶ موجب بهبود تخمین 
مقاومت فروروی خاک از طریق افزايش 1 از ۰/۷۵ به ۰/۸۰ و کاهش 
۳ به میزان ۸/۶۵ درصد گردید. سانتوس و همکاران ( 520105 
1 ,۸ ۶) بیان کردند که مقاومت فروروی خاک با دقت بالایی 
(۰/9۸ -*1 و ۰/۰۸۴ -808517) توسط شبکه عصبی مصنوعی از روی 
متغیرهای ورودی جرم مخصوص ظاهری و رطوبت خاک مزرعه مورد 


زین 


050 


061 


014 


اصغری و همکاران. مدل‌سازی مقاومت فروروی خاک با استفاده از رگرسیون. شبکه عصبی مصنوعی و برنامه‌ربزی بیان ژن 


تخمین قرار گرفت. شیخ‌زاده و همکاران ( ,اه ۶ 12۵060(ع٩‏ 
6 نیز در خاک‌های زراعی دشت اردبیل» مدل با متغیرهای ورودی 
رطوبت خاک مزرعه» سیلت» جرم مخصوص ظاهری و انحراف استاندارد 
هندسی (ع0) ذرات خاک را به‌عنوان بهترین مدل ۵۳۲ (52 ۲5۳ 
و ۷1۳ به‌ترتیب برابر با ۰۰/۴۱ ۳2 ۰۱/۴ 1۷۲۳۵ ۰/۲۰۶ - در داده‌های 
آزمونی) برای برآورد 18 خاک معرفی کردند. 


نتایج مدل‌های برنامه‌ریزی بیان ژن (۸(۷۲) برای برآورد 
مقاومت فروروی (۳18) خاک 

در برآورد ۳12 با استفاده برنامه‌ریزی بیان ژن (2۳)) نیز ۳ نوع 
مدل با همان متفیرهای ورودی مدل‌های ۷1 و ۸۱۷۱ (جداول ۵ و 
۶) ایحاد گردید؛ جدول ۷ معیارهای ارزیابی ۲ مدل 0۳۳ ارائه شده 
برای برآورد ۳1۶ را نشان می‌دهد. رتبه‌بندی مدل‌های ۳) نیز همانند 
مدل‌های ۸1,1 (جدول ۵) و ۸۳ (جدول ۶) براساس پایین بودن 
15 و بالا بودن ۸۷5 انجام گردید. بر مبنای آماره‌های مذکور در 
داده‌های آزمونی (جدول 4۷ دقیق‌ترین مدل 07۳ برای برآورد ۳15 
مدل شماره ۴ با متغیرهای ورودی رطوبت خاک مزرعه» سیلت و جرم 


۳۷۹ 


مخصوص ظاهری نسبی در مقایسه با سار مدل‌های 17۳) می‌باشد. 
مقادیر 422 215۳ 245 و 5 برای این مدل براساس داده‌های 
آموزشی به‌ترتیب برابر ۰/۷۵ ۷۳۵ ۰/۴۳ ۳۵ ۰/۰۴ و ۰/۷۴ و داده 
های آزمونی برابر ۰۰/۷۹ ۷۳۵ ۰۰/۹۱ ۷۳۵ ۰/۱۳ و ۰/۶۳ به‌دست آمد 
(جدول ۷). مقادیر مثبت 1۷11 در داده‌های آموزشی و آزمونی نشان‌دهنده 
اتکی مشص ورد کین تفای سا میل قو) شار:۲ 
می‌باشد. در رتبه‌های بعدی, به‌ترتیب مدل 1۳ شماره ۲ و شماره ۱ 
دارای بیشترین دقت در برآورد ۴ می‌باشد؛ با توجه به جدول ۷ مشاهده 
می‌شود افزودن متغیر سیلت به مدل 0۳۳ شماره ۱ باعث کاهش 
۳۴ از ۱/۲۶ به ۰/۹۳ مگاپاسکال و افزايش ضریب 38 از ۰/۲۸ 
به ۰/۶۱ در داده‌های آزمونی گردید. دیناروند و همکاران ( ۳027۷۵00 
8 ,۵ 2۶) نیز بیان کردند مدل 02 با متفیرهای ورودی رطوبت 
خاک مزرعه, جرم مخصوص ظاهری» رس و میانگین هندسی قطر (یه) 
ذرات خاک توانست ۳ خاک در اراضی زراعی دشت اردبیل با دقت 
بالایی تخمین بزند به‌طوری که مقادیر 122 ۷1517 و ۷170 مدل مذکور 
به ترتیب برابر با ۰۰/۵۱ ۳1۳۵ ۱ و ۷۲۳۵ ۰/۴۳- در داده‌های آزمونی 


به‌دست آمد. 


جدول 71- معیارهای ارزیابی مدل‌های شبکه عصبی مصنوعی مقاومت فروروی (۳1۸) خاک 
۴ 8 ۱۵06۱6 ۸ ۲۵۲ فان طمتاهت ۲۳۷۵ - 6 12016 


آموزش آزمون 
1۳2 ]1 16511۳02 
۹ 1۷1۳ ۱/۱ ۳ 5 11۳ ۳ ۳ متغیرهای ورودی شماره 
(۷۲۳۵) (۷۲۳۵) (۷۲۳۵) (۷۳۵) تا ۵ ۱۱۱ م۱( 
013 0.07 0.9 ۵4 0.43 017 112 06 1۳۳۷ 1 
0.70 0.3- 06 ۰0 0.83 09- 061 09 له ,۳۱۷۷ 2 
0.18 005 039 9 <۵ 092 5- 041 092 ۳1 مره ,۳۲۷۷۵ 3 


3 ضریب نش-ساتکلیف, 1۷8۳7: مجذور ميانگین مربعات خطاء :۷۲۳: میانگین خطاء؟: ضریب تبیین. توصیف متفیرها در زیر جدول ۲ آورده شده است. 
وعاصمننه ۷ عناوم منک مها :8 زعمجه صفعهد :۷۲ بتمجه مممنوه صقعصص نموم ۳6۴۰ +صمتعمتصمعاعل ۵۶ اصعئم‌تلاهمع :۳2 
۰ ۲20186 06۲صنا ماطاهاته2۷ ود ممتامتهوع 


جدول ۷- معیارهای ارزیابی مدل‌های برنامه‌ریزی بیان زن مقاومت فروروی (۳18) خاک 
۴ 0۶ ۱۵06۱ 6۲ 0۲ وتهانی صمتامتاه1۳7 -7 12016 
آزمون 


آموزش 
۵ ۱( 
11۳ ,9۳ 2 
(0۷۳۵) (۷۳۵) ی ۹ 
002- 059 05 028 
3- 053 061 061 
004 043 015 ۳ 


۷ 2و6 ] 
1۷۳ ,۳ 2 متغیرهای ورودی شماره 
(۷۳۵) (0۷۳) تا ۱۱۱۵۹۱۵۵۵۱ ۱/3 
008 126 030 ۵ 1 
013 093 ۳2( 1 ,۳۱۷۷۵ 2 
013 091 079 امرساظ لته ,۲۷۷ 3 


98لا ضریب نش-ساتکلیف, ۷5[7[؟: مجذور میانگین مربعات خطاء (۷1۳: میانگین خطاء 27 ضریب تبیین. توصیف متفیرها در زیر جدول ۲ آورده شده است. 
عناوم زامن مه :6 زرمجه صفعهد :۲۴ بتمجه متمنیوه صقعصص اممع ۳6۴۰ +صمتعصتصماعل ۵۶ )صوئم‌تزاهمع :۳2 
۰ ۲12016 06۲جنا ماصمأنته27 وا صمتاورتن‌ععل ععاماهته ۲۷ 


۰ نشریه آب و خاک. جلد ۳۸. شماره ۲. خرداد - تیر ۱۴۰۳ 


اندازه گیری۷]6۵50۲60 سره 
برآورد0002)60] صع حه 


مقاومت فروروی 
(0۵ 2۳ 
دبا 


۲1 ۲2 ۳3 ۳4 ۳5 ۳6 ۶1 ۶2 83 84 185 186 01 02 03 ۵4 5 06 07 8 


نقاط آزمون 


عاصلمم وصتاوع ]1 
شکل ۳- بررسی عملکرد بهترین مدل برنامه‌ریزی بیان ژن (612۳) در برآورد مقاومت فروروی (1) خاک در نقاط آزمونی (20-)؛ :: کاربری 
جنکلیء *1: کاربری مرتعی و ): کاربری زراعی 
عاصم عصتاعع) ما 2 (خ۳) معصه)عایی ممتاصامصهن (زمد ماهتا (م0م ای ۵ معصمممرمزهم عمط ممتاععتامی۲۸۲] -3 ۲12۱۷۲۵ 
۶6 ۱20 0۱0۱۲۵۲۵0 :2 2۳00 عون ۱2۳0 ۲۵2۳۵۵ :1 وکا 12۳00 10۲6۵8۲ ۲۰ :(20) 


بررسی عملکرد بهترین مدل در برآورد مقاومت فروروی (۳18) خاک 
در نقاط آزمون (۲۰ داده تصادفی) به‌عنوان نمونه در شکل ۳ برای مدل 
برنامه‌ریزی بیان ژن (012۳)) نشان داده شده است؛ با توجه به این شکل 
مشاهده می‌شود خطای مدل (اختلاف بین مقادیر اندازه‌گیری شده با 
مقادیر برآورد شده) در پیش‌بینی متغیر وابسته (۳1) از روی متفیرهای 
مستقل (رطوبت خاک مزرعه. سیلت و جرم مخصوص ظاهری نسبی) 
در نمونه‌های خاک برداشته شده از کاربری مرتعی (15) به مراتب بیشتر 
از کاربری‌های جنگلی (*) و زراعی (6) می‌باشد؛ علت این مسأله را 
شاید بتوان به ضریب تغییرات (۷)) بالای مقاومت فروروی خاک در 
کاربری مرتعی (۵۸ درصد) در مقایسه با کاربری‌های جنگلی (۲۲ درصد) 
و زراعی (۲۳ درصد) نسبت داد همچنین از بین متغیرهای مستقل نیز 
رطوبت خاک مزرعه دارای بیشترین ضریب تغییرات در کاربری مرتعی 
(۳۷ درصد) در مقایسه با کربری‌های جنگلی (۱۶ درصد) و زراعی (۱۶ 
درصد) بود. نتایج تحقیقات قبلی نیز بیانگر تأثیرپذیری بیشتر مقاومت 
فروروی از رطوبت خاک مزرعه در مقایسه با سایر ویژگی‌های زودیافت 
خاک می‌باشد ( 520605 :2011 ,.۵ 21 ۷2 :2008 .۵1 6 ]2ظ 


تتحقطوعه +2016 ره ۶ طمله2 معط :2011 ره 1 
9 ,0020201). 


مقایسه دقت مدل‌های 1۷۲۲ ۸ و 17۳ در برآورد ۳1٩‏ 
مقادیر 82 15۳ و 35 برای مدل‌های 11 (جدول ۵ 
۸ (جدول ع) و 02۳ (جدول ۷) ارائه شده برای برآورد ۳8 در 
داده‌های. آزموتی. پیانگر آن. است. که مدل‌های, شبکه‌های, عضبین 
مصنوعی (۸(۱) به علت داشتن 82 بالا (از ۰/۵۶ تا ۰/۹۲) و کل 
بیشتر (از ۰/۴۳ تا ۰/۹۳) و 1۳1510 پایین (از ۰/۴۱ تا ۱/۱۲ مگاپاسکال) 


در مقایسه با مدل‌های رگرسیون خطی چندگانه (7 از ۰/۲۷ تا ۰۰/۳۴ 
5 از ۰/۲۰ تا ۰/۳۶و 12/517 از ۱/۲۵ تا ۱/۳۲ مکاپاسکال) و مدل 
های برنامه‌ریزی بیان ژن (*12 از ۱۰/۳۰ ۰/۷۹ "٩‏ از ۱۰/۲۸ ۰۱۶۲و 

از ۰/۹۱ تا ۱/۲۶ مگاپاسکال) دارای دقت بالاتری و خطای 
کمتری در برآورد ۳ بودند. شکل ۴ نیز نمودار ۱:۱ مقادیر اندازه‌گیری 
و برآورد شده *۳1 خاک را با استفاده از بهترین مدل ,11 (مدل شماره 
۲ در جدول ۵ بهترین مدل ۸ (مدل شماره ۳ در جدول ۶) و 
بهترین مدل 012۳ (مدل شماره ۳ در جدول ۷) برای داده‌های آزمونی 
نشان می‌دهد. با توجه به این شکل مشاهده می‌شود که پراکندگی داده 
ها حول محور ۱:۱ در بهترین مدل ۸۷۱ کمتر از مدل‌های ۷11 و 
لت برای برآورد ۳ در خاک‌های منطقه مورد مطالعه می‌باشد 
همچنین می‌توان گفت مدل‌های شبکه عصبی مصنوعی و برنامه‌ریزی 
بیان ژن در تخمین ۲15 به‌دلیل استفاده از روابط غیر خطی به مراتب 
دقیق‌تر از مدل‌های رگرسیونی خطی چندگانه عمل نموده‌اند علت 
برآورد دقیق‌تر شبکه‌های عصبی مصنوعی نسبت به مدل‌های 
رگرسیونی را می‌توان در نحوه عمل شبکه پرسپترون سه لایه دانست. 
عملکرد این شبکه به این صورت است که ابتدا الگویی به شبکه عصبی 
عرضه می‌شود و خروجی آن محاسبه می‌شود سپس شبکه این خروجی 
را با مقادیر اندازه‌گیری شده متغیر یا متفیرهای خروجی مقایسه می‌کند 
و ضرایب وزنی شبکه را به نوعی تغییر می‌دهد که کمترین اختلاف بین 
مقادیر پیش‌بینی شده و اندازه‌گیری شده باشد. بیات و همکاران ( ۲2726 
8 ,.۵1 2) و شیخ‌زاده و همکاران (2016 ,باه 6 1۳2۵06۲عط) 
نیز گزارش نمودند که شبکه‌های عصبی مصنوعی نسبت به مدل‌های 
رگرسیون خطی چندگانه در برآورد ۳18 دقیق‌تر عمل نمودند. سانتوس 
و همکاران (2011 ,.۵1 6 ۹20105) نشان دادند تابع شبکه عصبی 
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مصنوعی (۰/۹۸ 1 و ۰/۰۸۴ -81512) نسبت به تابع رگرسیونی های رگرسیون خطی و غیرخطی در برآورد ویژگی‌های دیریافت خاک 
(۰/۹۲ 2 و ۰/۹۵۱ -۳(۷/۹۳) دقت بیش‌تری در برآورد مقاومت از روی ویژگی‌های زودیافت خاک توسط احمدزاده کلیبر و فولاپناه 
فروروی از روی جرم مخصوص ظاهری و رطوبت خاک مزرعه داشت. (2023 مطممهم‌ن ۳۵۱20‏ تدطاتم‌ل۴۵ ط02۵0مصطض) و زانگ و 
نتایج مشابهی در خصوص دقت بالای مدل طتآن) در مقایسه با مدل زانگ (2024 ,طععط2 ک عحمطم2) گزارش گردیده است. 
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(۲۳۵) ۳ ۷۲625۱۲60[ 
مقاومت فروروی اندازه گیری شده 


شکل 6- مقایسه مقادیر اندازه‌گیری و برآورد شده مقاومت فروروی (۳16) خاک با استفاده از بهترین مدل رگرسیون خطی چندگانه 
(۷]1» شبکه عصبی مصنوعی (((۸) و برنامه‌ریزی بیان زن (217۳)) بر اساس داده‌های آزمونی (20-). اطلاعات مربوط به مدل‌ها و متغیرهای 
ورودی در جداول ۰۵ ۱ و ۷ آورده ده است. خطچین بیانگر نمودار ۱:۱ است. 
6 جه 60عد )56 0262 هصتامع) 6 ۲۵۲ (۳۱۵) مصه)مامی۲ نموه ۵ کما۲۷۵ تامهم 0اه 60)مصصنای ۵ مصرمعاتوم‌صمن -4 عاع۲1 


6 ۲۵۲۰۱۵0۱۵6۹۰ )نامرد صرح عم جمهع )مصاه مت فص 062260 ۵۲۱۵6 6 20 5 ر4 م1201 .عا06مص نی 200 ۸۱۷ بخلن]۷ اععط 
۰ 1:1 ۲06 1016206 معصنا 0مصاعول 


۳ شریه آب و خاک. جلد ۳۸. شماره ۲. خرداد - تیر ۱۴۰۳ 


نتیجه گیری کمتری در برآورد ۴ خاک نسبت به مدل‌های ]۷۲ بودند. پیشنهاد 
می‌شود در تحقیقات بعدی از روش رگرسیون غیر خطی و نیز سایر 
مدل‌های هوشمند مانند نروفازی» ماشین بردار پشتیبان و جنگل 
تصادفی برای برآورد مقاومت فروروی خاک استفاده گردیده 9 نتایج آن 
با مدل‌های رگرسیونی مقایسه گردد. 


در این تحقیق از مدل‌های رگرسیونی خطی چندگانه (11 
شبکه عصبی مصنوعی ((۸(۷) و برنامه‌ریزی بیان ژن (طت1)) برای 
برآورد مقاومت فروروی (۳۷۹) خاک از روی ویژگی‌های زودیافت خاک 
در اراضی جنگلی» مرتعی و زراعی فندقلوی اردبیل استفاده گردید. نتایج 
1۳۳۳ 
۱ ِ و مس این مقاله مستخرج از طرح پژوهشی مصوب دانشگاه محقق اردبیلی 
ویژگی‌های زودیافت خاک از اهمیت بیشتری برخوردار بوده و در ورودی ۱ 
مدل‌ها ب‌کار گرفته شدند. مقایسه نتایج مدل‌ها نشان داد ۸0۱۷ پر می‌باشد به این وسیله از حمایت‌های معلونت محترم پژوهشی و فن 
علت داشتن ضرایب تبین (*8) و نش-ساتکلیف (۵15 بالات و میعنور ۰" آوری دانشگاه تشکر و قدردانی می‌شود 
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